Determining the accuracy of spatial gradient sensing using statistical mechanics 
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Many eukaryotic cells are able to sense chemical gradients by directly measuring spatial concen- 
tration differences. The precision of such gradient sensing is limited by fluctuations in the binding 
of diffusing particles to specific receptors on the cell surface. Here, we explore the physical limits 
of the spatial sensing mechanism by modeling the chemotactic cell as an Ising spin chain subject 
to a spatially varying field. This allows us to derive the maximum likelihood estimators of the gra- 
dient parameters as well as explicit expressions for their asymptotic uncertainties. The accuracy 
increases with the cell's size and our results demonstrate that this accuracy be further increased by 
introducing a non-zero cooperativity between neighboring receptors. Thus, consistent with recent 
experimental data, it is possible for small bacteria to perform spatial measurements of gradients. 
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Cells often direct their motion under the guidance of 
chemical gradients. This is essential for critical biological 
functions including neuronal development, wound repair 
and cancer spreading [HQ]. To detect gradients, small 
organisms like bacterial cells usually employ a tempo- 
ral sensing strategy by measuring and comparing con- 
centration signals over time along their swimming tracks 
In contrast, eukaryotic cells are sufficiently large 
to implement a spatial sensing mechanism, as they can 
measure the concentration differences across their cell 
bodies. Measurements for both strategies are accom- 
plished by specific cell-surface receptors which diffusing 
chemical particles (ligands) can bind to. Spatial sens- 
ing among eukaryotes exhibits a remarkable sensitivity 
to gradients of merely 1-2% across the cell . Given 
the dynamic fluctuations in ligand-receptor interaction, 
the receptor signal is inherently noisy, as demonstrated 
by single-cell imaging experiments [loL [ll[ . This natu- 
rally raises a question concerning the reliability of spatial 
gradient sensing. 

In analyzing bacterial chemotaxis, Berg and Purcell 
showed that the minimal uncertainty of mean concentra- 
tion measurements is set by the diffusion of ligand parti- 
cles [HJ ■ This work has been extended to include ligand- 
receptor binding effects and possible receptor coopera- 
tivity [l3|-[I3 • The spatial sensing program concerns the 
acquisition of information regarding the asymmetry in 
space (the gradient steepness and direction). The accu- 
racy of gradient measurements should similarly be lim- 
ited by physical laws governing diffusion and stochastic 
ligand-receptor dynamics. There have been some studies 
on the limits to spatial sensing, but either for idealized 
sensing mechanisms that ignore the receptor kinetics [l6[ 
or based on specific transduction models [HI . In this Let- 
ter, we address this problem in a general way using a sta- 
tistical mechanical approach where we view the surface 
receptors as a (possibly coupled) spin chain and treat the 
chemical gradient as a perturbation field. By calculat- 
ing the system's partition function, we are able to derive 
the physical limits of gradient sensing for both indepen- 
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FIG. 1: Schematic representation of our model: a circular cell, 
covered with receptors, is placed in an exponential gradient. 
The forward and backward rates k± control the transition 
between the bound and unbound states for the receptors. 



dent receptors and for receptors exhibiting cooperativity. 
These limits allows us to predict that the strategy of spa- 
tial sensing may not be exclusive to large eukaryotic cells 
but may also be applicable to some bacterial cells [19j . 
especially with the aid of receptor cooperativity. 

We consider a circular cell with diameter L immersed 
in a chemoattractant gradient (Fig. 1). We suppose that 
there are N receptors distributed at equally spaced in- 
tervals on the cell's perimeter [2(|. The angular coordi- 
nates of these receptors are indicated by ip n = 2ivn/N for 
n = 1, ...,N. For analytical convenience, we assume that 
the gradient field takes an exponential profile, as was re- 
cently realized in exp eriments utilizing the social amoeba 
Dictyostelium [9|, [21| . The local concentration at the n-th 
receptor is C n = Co exp [| cos(ip n — <f)\ , where Cq is the 
background concentration, p = (j-^cr denotes the gra- 
dient steepness, and <f> indicates tne gradient direction. 
Like a spin in physics, each receptor switches between 
two states: either active (s n = +1) or inactive (—1). For 
independent receptors, a receptor is activated only if it is 
bound by ligand and inactive otherwise. Let the energy 
associated with the state s n = +1 (or —1) be — e n (or 
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+e„) in units of the thermal energy /c^T. Then the "on" 
probability of the n-i/i spin is given by the Boltzmann dis- 
tribution: P on — e e "/(e £ " + e~ £n ). For simple receptor- 
ligand kinetics (Fig. 1), we have P on — C n /{C n + Kd) 
in chemical equilibrium where K d = k-/k + is the dis- 
sociation constant. Therefore, the free energy has the 
expression: 
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+ -cos(ip n -4>) = a +h n . (1) 



We define three statistical quantities (zq,Zi,Z2) = 
cos tp n , g J2n s ™ s ^ nl Pn) where z is a 
measure of the average receptor activity and where Z\ 
and Z2 measure the asymmetry in the receptor activ- 
ity. Using the transformation (01,02) = (p cos 4>, P sin c6) 
we can write the system's Hamiltonian as 'Hjv{sn} = 
-^ n £„s„ = — ao^o — (ctizi + a 2 z 2 )/2 and compute its 
logarithm partition function as follows, 
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where in the last step the summand is expanded in powers 
of p and the sum is replaced by an integral over [0, 2ir}. 

The partition function contains all the thermodynamic 
information we need to infer the gradient parameters p 
and </>, or alternatively, the transformed parameters a\ 
and a 2 . Since p 2 — a\ + a 2 ., we have by Eq. (2): 
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In addition, one can check that Cov[zi,Z2] = 0. Thus, 
for small p, the joint probability density of z\ and z 2 is 
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with fi = NC K d /(4(C + K d ) 2 ) and a 2 = 2fi \2$. It 
is easy to show that the maximum likelihood estimator 
(MLE) [22| of ai )2 is 812 = ^1,2//^- As an orthogonal 
transformation, the MLE of p and <p are given by p = 
\/a 2 + a 2 = \/ z\ + z 2 and <fi — arctan(a2/ai) = 
arctan(z 2 /zi), respectively. By the properties of MLE, 
both p and <fi tend to be unbiased and normal in the large 
N limit, i.e., p A N(p,cr 2 ) and 4> — * N(<f>, cr?), where 

"— denotes convergence in distribution. The asymp- 
totic variances in the gradient steepness and direction, 
' can be derived from the Fisher information 
which is diagonal since p and cj> are indepen- 
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FIG. 2: (CcA m OiklilKi): i(A) T- ftc irneq&amlie^ a 1 and ai 
versus ln(C /K d ) (X ( ;^'= 10%, ^ p = 8 WWoYand versus p 
(B; Co = K d , N = 80000). The solid lines correspond to 
the approximate analytical expressions while the symbols are 
the result of 5000 independent Monte-Carlo realizations. (C) 
<r| as a function of ln(Co/Kd) for different values of J. (D) 
The critical cell size below which spatial gradient sensing is 
ineffective, normalized by the critical cell size in the absence 
of cooperativity, as a function of the cooperativity strength. 
In C and D, we have chosen N = 80000 and p = 8%. 



dent. Thus, 
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and thus er 2 = cr 2 /p 2 . According to the Cramer- Rao 
inequality, a 2 and cr? set the lowest uncertainties of gra- 
dient measurements from an instantaneous sampling of 
the receptor states [23 |. The approximation for both 
variances is plotted in Fig. 2 as a function of the two 
parameters characterizing the gradient: the background 
concentration Cq (Fig. 2A) and the gradient steepness p 
(Fig. 2B). We have also performed Monte-Carlo simula- 
tions in which 80000 receptors are uniformly distributed 
along the circular cell membrane. Computing our statis- 
tical quantities for 5000 independent realizations, we de- 
termined the measurement errors and have plotted them 
as symbols in Fig. 2. The analytical results agree well 
with the numerically obtained values. From Fig. 2, we 
can see that the variances reach a minimum for Co = K d 
while only the error in the gradient direction depends on 
the steepness of the gradient (er 2 , ~ p~ 2 ). Thus, since 

p = pqL (with po — (j-^r) increases with the cell's size, 
larger cells are able to sense the gradient direction with 
higher accuracy. 

The above results are derived from a single snapshot 
of the system. If the cell integrates receptor signals over 
some time interval T, then averaging over multiple mea- 
surements can appreciably reduce the errors of gradient 
sensing. However, the capacity of such averaging is lim- 
ited by the expected time it takes for every independent 
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measurement. As shown in [§1, the time to com- 
plete a single measurement is roughly twice the system's 
correlation time r resulting from the diffusion and bind- 
ing of ligand molecules, leading to a reduction of the 
variance: &pj- — ^f°p- The correlation time is given 
by t — r rcc + T diff , where r rcc = l/(fc_ + C k + ) is the 
timescale of receptor-ligand reaction and Tdiff describes 
the diffusive transport time of ligands. Let 77 = Tdiff /T rcc , 
then the measurement is reaction-limited if rj <C 1 and 
diffusion- limited if r\ 3> 1. From the above arguments we 
find that averaging signals over T yields a lower uncer- 
tainty of the gradient estimate, 
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For small background concentrations (Co <C Kd), Tdiff = 
N / {2irLDKd) where D denotes the ligand diffusion coef- 
ficient [HI, Q, , and the uncertainty reduces to a 2 j- ~ 
16/ (NTCok + ) + 8/(nTDLCo). This expression contains 
two terms: the first one is determined by the chemical ki- 
netics, and the second one, up to a geometric constant, is 
exactly the Berg-Purcell limit [12j or the result recently 
derived in [ll| . We can derive similar results for the direc- 
tion inference, since u 2 j- = cr 2 t/p 2 - For typical eukary- 
otic cells, it has been estimated [1J,[23| that 77 <C 1, which 
implies a\ r ~ 16(1 + K d /C )/(Np 2 Tk.) ~ l/(Np 2 ). 
We can relate the uncertainty in the direction measure- 
ment to the cell's size L. Assuming that the number of re- 
ceptors in our model scales with the cell size as N = NqL s 
with ^ 6 ^ 2, we find a 2 ^ j- ~ L~( 2+A ). For comparison, 
the Berg-Purcell analysis considered only an average con- 
centration measurement and scales as o c j- ~ 
Not surprisingly, our results indicate that spatial direc- 
tional sensing can be more sensitive to the cell's size. 

Our analysis above, which extends beyond the Berg- 
Purcell framework by providing a direct calculation of 
the directional sensing limit er^ 7-, was carried out for 
independent receptors, as is assumed to be the case for 
most eukaryotic cells that have been studied to date. We 
now ask, what if there is receptor coo per ativity as has 
been found in many bacterial cells |2J-[2a)? Intuitively, 
short-range interactions make it possible for receptors 
to collectively respond and thus sharpens the asymme- 
try of receptor signals. It is natural to speculate that 
such enhanced sensitivity may set new and lower limits 
for directional sensing. To incorporate potential receptor 
cooperativity, we extend our model to include a nearest- 
neighbor interaction J (again, in units of the thermal 
energy ksT). Now, the activity of a receptor, again rep- 
resented by {s n = ±1}^ =1 , is determined not only by 
the local chemical concentration but also by the states of 
its neighboring receptors. This means that an unbound 
receptor is not necessarily inactive, as it may have been 
affected by active, nearby receptors. 

Because the local concentration is identical for nearest- 
neighbor sites (i.e., e n — e n ±i), the Hamiltonian of 
our Ising chain can be written in a symmetric form: 

H. N {s n } = -Elll^A+l + £n{Sn + S n+1 )/2], with 



the boundary condition sn+i = si- The correspond- 
ing partition function is Qn = J2 S —Y^ SN e~( Ho+Hl >, 
where H = -J2 n \- Js n s n+l + a (s n + s n+1 )/2] repre- 
sents the Hamiltonian of an isotropic reference system 
and where Hi = -J2n s n h ™ = "4 En 5 "^" _ 4>) 
results from the spatial heterogeneity of the concentra- 
tion. For small p, one can view Hi as a perturbation to 
H . The partition function of the reference system, Q^' , 
is exactly solvable (27j, e.g., using the transfer matrix 
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with A± = e J coshao ± y e 2J + e 2J sinh 2 ao being the 
eigenvalues of V . Thus, In Q^} — > N In A+ for large N, 
The statistical perturbation theory inspires us to write 



Qn - Q^n J2si - ^s 



-H 



Hl IQ { n = Q% } {e~^) 



+ 4 En( S ") C0S6»„ -I- §2 I]„, m (SnS m ) COS 6»„ COS m ]. 

Here, we denote 6 n = tp n — 4> for short and use (•) to rep- 
resent the expectation over the reference system. Due 
to isotropy, (s n ) is independent of its location (index n) 
and hence J2n( s n) cos9 n = (s n ) Y^ n cos ^« = 0. We fur- 
ther calculate that ^2 n m (s n s m ) cos 9 n cos 9 m — ^(1 + 

20/(1 + e 4J sinh 2 a ), fi], where £ = [ln^/A.)]" 1 de- 
fines the correlation length of the classic Ising chain [13] ■ 
Finally, the log-partition function of our model is 
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which reduces to Eq. (2) as J — > 0. 

Now we rewrite Hn = —Jj2 n s nS n +i — ao-zo — {a\zi + 
&2Z2)/2, with the same notations for a% and Zi, i = 0, 1, 2. 
As has been demonstrated before, the MLE of ai and a-i 
can be found from the joint Gaussian distribution of zi 
and Z2, except now we have to replace /ibypt = j^N(l + 
2£)/(l+e 4J sinh 2 a Q ). So the MLE oip and (f> are given by 

p = 7i _1 y< ' z\ + z\ Aa/"(p,?p) and 4> = arctan(z 2 /^i) — > 
Af((f>, ct^). Similar to Eq. (5-6), their variances are a 2 = 
2//I and a\ = a 2 /p 2 = 2/(Jlp 2 ) We plot a\ as a 

function of \n(Co/Kd) for different values of J in Fig. 
2C. Regardless of the receptor coupling strength, this 
error is minimal at Co = K d (or ao = 0) where the 
correlation length is £ = l/ln(coth J) ~ \e 2J and ct 2 , ~ 

32/[Ap 2 (l + e 2J )] -^/(l + e 2J ). 

Receptor cooperativity may help a smaller cell of diam- 
eter L achieve the same level of accuracy as a larger cell 
of diameter L with independent receptors, i.e., ?^(L) = 
a 2 (L). By our previous scaling assumption, the recep- 
tor number of the smaller cell is N = N(L/L) S . If L* 
denotes a critical cell length below which spatial sensing 
is infeasible with non-cooperative receptors, then adding 
cooperativity will push the critical cell size, L* , lower by 
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a factor of (1 + 2f)- 1/(2+5) ~ L(l + e 2-/)-i/(2+5). This 
is shown in Fig. 2D where we have plotted L*/L* as 
a function of J for three values of the scaling factor i5. 
As a specific example, we take L* — 8fim which cor- 
responds to the typical size of a Dictyostelium amoeba. 
Then, we see that for a cooperativity of J = 0.5 the new 
critical size becomes L* ~ 4 — 6/im, comparable to the 
size of many bacterial cells. It is worth remarking that 
although receptor interaction improves the precision of 
gradient sensing for Cq close to Kd, it enlarges the er- 
rors when Co is far away from Kg (Fig. 2C). In other 
words, the improved accuracy near is at the cost of 
the sensitivity range of background concentrations. Such 
a tradeoff could be a limiting factor for the introduction 
of coupling into the spatial sensing mechanism. 



It is commonly believed that prokaryotic cells such as 
E. coli are too small to perform spatial sensing of chemi- 
cal gradients. However, recent experimental observations 
show that at least one type of vibrioid bacteria (typical 
size 2x6 /im) are able to spatially sense gradients along 
distances as short as 5 [ira [19| . Our results allow for the 
possibility that smaller organisms employ a spatial sens- 
ing strategy with the aid of receptor cooperativity. As 
spatial sensing is argued to be superior to temporal sens- 
ing for fast swimming bacteria [3(| , this possibility is 
of significant theoretical interest and remains a challenge 
for future empirical studies. 
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